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STUDIES 


ON  THE  SEPARATION  OF  FLOW  AT  THE  EDOES  OF 
ANTI SYMMETRICAL  DELTA  WINGS  IN  SUPERSONIC 
SYSTEM 

by  Elie  Carafoll  and  Stefan  Staicu 

Institute  of  the  mechanics  of  fluids 

at  the  Academy  of  the  Socialist  Republic  of  Romania 

In  this  work  we  will  study  the  supersonic  flow  around 

anti symmetrical  delta  wings,  taking  Into  consideration 

the  separation  of  flow  at  the  subsonic  leading  edges. 

As  with  the  thin  delta  wing  with  constant  Incidence,  In 

this  case  the  flow  separates  around  the  edges  creating 

a vortex  layer  which  is  situated  on  the  upper  side  of  the 

wing  transforming  Itself  Into  two 

of  the  same  Intensity  and  sign,  having  however  q-t*- 

anti symmetrical  position  the  axis  of  symmetry  of 

the  wing.  The  formation  of  vortexes  directly  effect  the 
IcomplexJ 

production  of  avf laid  of  vertical  velocity,  which  will 
modify  the  flow  In  such  a way  that  the  pressures  will 
be  finite  at  the  leading  edges.  The  distribution  of 
vertical  velocities  leads  to  a system  of  three  Imaginary 
wings,  through  the  superposltlonlng  of  which  a resulting 


ivortexi 

cone  en  t r a tedrnucl  eu  s e s 
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imaginary  wing  Is  obtained,  equivalent  from  an  aerodynamic 
point  of  view  with  the  real  wing.  The  flow  which  occurs 
is  homogenous  from  the  second  order,  and  we  have  the 
possibility  of  determining  the  pressures  on  the  upper 
and  lower  sides  of  the  wing  and  the  aerodynamic  characteristics. 

1.  Preliminary  Considerations 

In  the  present  work  we  will  study  the  supersonic  flow 
around  anti  symmetrical  thin  delta  wings,  taking  into  consideration 
the  falling  off  of  flow  around  the  subsonic  leading  edges. 

The  antisymmetrlcal  distribution  of  the  incidences  or  of  the 
vertloal  velocities  corresponds  to  a torsioned  delta  wing 
antisymmetrlcal  according  to  a llniar  function.  The  velocity 
of  the  undisturbed  flow  to  be  Uoo  , parallel  with  the  axis  of 
symmetry  of  the  wing  (Fig.  1). 

As  with  the  thin  delta  wing  plane  with  constant  incidence 
(1),  the  flow  separates  at  the  leading  edge,  creating  a vortex 
layer  situated  above  as  well  as  below  the  wing,  producing  an 
antisymmetrlcal  movement.  The  vortex  layer  having  sufficiently 
small  thickness  can  be  considered  like  a vortex  sheet,  which 
winds  Itself,  as  various  authors  have  shown,  among  whom  we 
will  mention  M.  Roy  (4),  in  the  form  of  a horn  composed  of  a 
concentrated  nucleus  and  a marginal  vortex  sheet  starting  from 


J 
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the  leading  edge 


Since  the  Incidence  Is  variable  on  the  surface  of 
the  wing,  the  axis  on  which  the  horn  winds  will  be  a curve 
and  the  vortex  generation  Intensity  of  the  nucleus  is 
variable  along  the  axis,  proportional  with  the  span  of  the 
wing.  For  simplification.  In  the  following,  the  axis  on 
which  the  nucleus  of  the  vortex  Is  situated  Is  considered 
to  be  a straight  line. 

In  this  manner  the  field  of  flow  Is  modified  by  the 

nucleuses 

existence  of  two  of  vortexes  of  the  same  Intensity 

and  sign,  situated  antisymmetrically  In  Reference  with  the 
axis  of  symmetry  Oxj  (Fig.  l)  at  the  abscissa  c and  the 
ordinate  t. 


Making  these  considerations,  the  flow  around  the  wing 
remains  In  continuation  conical  of  the  second  order  and  can 
be  treated  using  methods  from  the  theory  of  conical  motion 
of  the  higher  order  (2).  Thus,  for  the  solution  of  this 
problem,  we  will  follow  the  way  used  In  previous  papers  (1), 
(3)t  (11)»  where  solutions  were  given  for  the  thin  delta  wings, 
with  constant  and  anti symmetrical  Incidence,  respectively, 
(forced  antisymmetry)  In  reference  to  Oxi,  which  led  us  to 
a conical  motion  as  a matter  of  fact. 


3. 


For  this,  we  will  allow  that  the  effect  of  the  falling 
off  of  flow  at  the  leading  edges  and  the  formation  of  the 
two  antlsymmetrlcal  nuclei  is  to  creat  a vertical  and  long- 
itudinal field  of  velocity  which  will  bring  about  modifications 
on  the  field  of  flow,  hmvlng  as  a result  the  avoidance  of 
infinite  velocities  at  the  leading  edges,  as  results  from 
classical  linear  theory.  But  it  can  be  allowed  that  the 
effect  of  the  longitudinal  velocities  of  disturbance  will  be 
able  to  be  substituted  through  that  of  a distribution 
corresponding  to  the  vertical  velocities.  For  this,  we  will 
consider  a distribution  of  vertical  velocity,  incidence  and 
antisymmetry  respectively,  so  as  to  correspond  to  a real  case 
of  an  imaginary  thin  delta  wing  with  variable  incidence, 
different  on  its  two  sides,  having  at  the  same  time  finite 
velocities  at  the  leading  edges.  Due  to  antisymmetry,  the 
additional  axis  and  vertical  velocities  created  by  the  nucleus 
of  the  vortexes  expire  toward  the  middle  of  wings,  becoming 
equal  with  zero  on  the  axis  of  symmetry  Oxj . In  this  way  the 
effeots  of  the  vortexes  on  the  central  line  are  canceled, 
where  later  a nought  vertical  velocity  will  remain.  We  will 
allow  that  the  real  thin  wing,  which  has  in  a certain  way 
finite  velocities  at  the  edges  through  the  effect  of  the  falling 
off  of  flow,  is  equivalent  from  an  aerodynamic  point  of  view 
with  an  imaginary  thin  wing,  having  the  same  variation  of 
incidences  which  we  defined  above. 
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In  order  to  study  the  motion  efcsily,  through  the  conical 
notion  method,  we  will  split  the  Imaginary  wing  corresponding 


to  the  distribution  of  vertical  velocity  above  Into  three 
wing  components,  as  we  have  proceeded  in  previous  cases  ( ll , 
(3)* 


1}  The  thin  wing  having  the  variation  of  anti  symmetrical 
Incidence  suitably  chosen  to  respect  somewhat  the  phenomena 
of  pressure  modifications  and  vertical  velocities  on  the 
surface  of  the  wing  near  the  leading  edge*  In  this  way  we 
obtain  an  Imaginary  thin  wing  with  finite  velocity  at  the  leading 
edges,  of  mqual  and  opposed  sign  on  the  two  faces,  higher  and 
lower. 


i 


2)  The  wing  of  symmetrical  "thickness"  - the  notion  of 
"thickness"  here  has  an  Imaginary  sense  - having  the  slope 

equal  and  of  the  same  sign  as  the  Incidence  of  the  first  wing. 
This  wing  combined  with  1)  will  have  different  pressures  on 
the  two  faces,  as  happens  In  reality. 

3)  The  third  wing  will  have  symmetrical  "thickness"  with 
variable  and  antlsymmetrlcal  slope,  however  such  that,  combined 
with  wing  2),  to  obtain  a mean  nought  slope,  corresponds  to 

a real  thin  wing.  Through  superposltlonlng  the  three  wing 
components,  we  will  obtain  the  resulting  imaginary  wing 
equivalent  from  an  aerodynamic  point  of  view  with  the  delta 
wlftg  with  the  separation  of  flow. 
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Examining  the  real  phenomena  of  the  falling  off  of  flow 
at  the  edges  taking  Into  consideration  experimental  Indications, 
we  will  manage  to  chose  a distribution  corresponding  to  the 
vertical  velocities  or  the  Incidences. 


I ‘j 


Tn  this  way  we  will  further  mark 

w':  - «:  t m,  . _ . p 

*’  (1) 

the  vertloal  velocities  and  the  Incidences,  w^,<x'u  , on  the 
higher  side  respectively,  on  the  lower  of  the  imaginary 

thin  wing}  vertical  velocity  » -■  Vm  on  the  real  wing  ia  given 

by  the  relation 
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2.  THE  DETERMINATION  OF  THE  AXIS  OF  DISTURBANCE 
VELOCITIES 


In  continuation  we  will  determine  for  the  three 
anti  symmetrical  wing  components,  the  axis  of  disturbance 
velocities,  being  necessary  for  the  calculation  of  the 
distribution  of  pressures  and  of  the  aerodynamic  characteristics 
of  the  resulting  Imaginary  wings,  which  are  presupposed  to  be 
the  same  as  real  thin  delta  wings,  having  the  Incidence 
defined  by  (2).  The  motion  around  the  wings  being  conical 
of  the  second  order,  we  will  again  use  the  methods  used  before, 
considering  In  this  sense  the  physical  plane  Oyz  (fig.  1) 
normal  on  Ox^  and  having  the  coordinates 


the  axises  Oy  and  Oz  being  parallel,  respectively  with  0x2 
and  Ox^.  Again  we  will  make  the  transformation  similar  with 
that  given  by  Busemann  (fig.  2)* 


•ii  = 




l - H*z* » 


S 


*Jff  - BHy*  -f  *») 
i - J 


(x  = D + ii), 


obtaining  a nlane  that  has  the  property  of  keeping  the  track 
of  the  wing  (*y»  y,  z ■ 0)  in  the  true  magnitude. 


As  we  know,  In  this  plan  the  first  derivatives  of  the 


velocities  of  disturbance  u,v,  w are  harmonic  functions  and 
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the  respective  conjugated  functions  can  be  associated  In 
such  a nanner  that  analytic  functions  of  coaplex  variables 
x ■ cJ  ♦ 1^  (4)  are  obtained* 


1)  the  thin  ant lsymmetrlcal  wing.  We  will  consider  the 
vertical  velocity  on  this  wing,  as  a result  of  the  effect  of 
the  two  vortex  nuclei,  constant  along  axis  Oxj  and  variable 
on  the  track  of  the  wing  according  to  the  relation 

«■'  = K'ui  (y)-r,  = -r,y  MV,  (y),  !r\ 


Ai 


t — r 

r 5 - 
r 7 


- f/fl 1 


In  which  the  parameter  WQ^(y),  variable  on  the  wing,  Is 
obtained  with  the  help  of  a distribution  corresponding  by 
singularity  of  the  second  order  (source)  placed  on  the  track 
of  the  wing  feoa  the  plane  x + 1^,  with  the  Intensity 
<{  given  by  the  function 

q'(y)  ==  «*(  1 - y)  (*  < y < *>,  /w 

we  obtain  In  this  way  the  velocities  -r,  -=  - r.x,/,  t at  the 
edge  of  wing  **  - *i1  ■ r.  .r,  y for  the  points  on  the  wing 
contained  between  y ■ o and  y ■ a* 


In  this  manner.  In  the  case  of  nomogenous  motion  by 
order  two  (n  - 2),  we  will  have  the  following  distributions 

of  the  Intensities  of  sources: 

vlo  1 1 - -y)'  ± «”(  1 - T")  ('  < T<  ^ l)’  (7> 

where  n will  represent  the  position  of  an  element  of 

intensity  on  the  track  of  the  wing  from  the  auxiliary  plane 

x „ 3 ♦ Keeping  In  mind  previous  works  (2),  the  contribution 

of  elementary  sources  dq'  In  the  expression  of  the  axis  of 

disturbance  velocity  from  point  x will  be 

,,v  , ...  I 0 4-n)U- jj.  (8) 

All  ^ | •Ji(T)—  j) 

The  axis  of  disturbance  velocity  for  the  thin  wing  component 
will  be  obtained  through  the  addition  of  all  the  contributions 
of  sources  under  the  form 


i * ^ ♦ HO  - 

-tK1  - 1)'*-  **  w 


which,  after  completing  the  calculations,  becomes 

t * — *x 


I A 

yjj- 


1 , l I*  — *X 

|i£  , _ -A-  f qu  J)(»  - J)  [ 1 ~ -yf  P*  h ((,  _ x) 


l * -x\  I I*  4 •*  , 

**,  - «/*,  *)  (*  + ')  ( 1 - ~yr)  °OS  ' |(«  1 .r) 

^ - 2?I«f)  ' -y  + fl!*1  - J.  ) \ 1 /*  1' 


(10) 
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2)  The  wing  of  symmetrical  thickness,  with  the  slope 
equal  with  the  Incidence  of  the  first  thin  wing.  The 
Introduction  of  the  second  wing,  of  symmetrical  "thickness". 
Is  necessary  In  order  to  remove  accentuated  peaks  of  pressure 
on  the  lower  side  of  the  wing,  where  the  division  of  the 
pressures  obtained  through  superpositioning  with  the  first 
wing  component,  will  be  different  from  that  on  the  higher 
side  of  the  wing.  Following  the  general  method  of  conical 
motions  (2),  for  a wing  of  symmetrical  thickness  with  the 
variation  of  slope  given  by  the  same  distribution  of  sources 
(7),  we  will  write  for  the  axis  of  disturbance  velocity  the 
following  expression,  similar  with  (9): 


-Hl'-i)1*-* ■(/" 


Completing  the  above  Integrals,  we  find 


■j-  — (Uit  - — 4 q*, 

/ \l  I * 4-  *1  . _ I 1 — ] 

l 21  ) HI*  - x)  J 


1 ~ 

m - x) 


- (,t  - + *>(•-  V)  iu?- 


10 


/ H — X\  . . 1 + ^**£.1 

- (•  + x)  ( 1 “ 2 1 ) C°8  1 + *)  1 

xV 1 - &***[  (Bill-*  Bl  - Bin'1  #*)  (</’«  - 2 ^ ~ 

iii  l 

— — (^l  — BW  — Vl  - H4*2)]}  ± 11 

where  L represents  the  contribution  of  the  subsonic  edge 
having  the  slope  equal  with  *®>l- 

3)  The  wing  of  symmetrical  thickness  compensating  for 
slope.  In  continuation  we  will  remark  that  through  the 
Introduction  of  the  effect  of  the  wing  from  point  2),  the 
resulting  wing  became  the  "large"  wing.  We  will  be  able  to 
compensate  this  work.  Introducing  on  the  track  of  the  wing 
from  plane  x (4)  a new  distribution  of  source  of  an  adequate 
form,  which  will  restore  the  wing  to  an  average  nought  thickness. 
The  variation  of  the  vertical  velocities  w'J'  given  by  these  sources, 
which  will  necessarily  present  at  the  edge  of  the  wing  the 
velocity  — will  correspond  with  a "wing  compensating  for 

slope"  of  symmetrical  "thickness",  which  has  the  role  of 
canceling  the  mean  slope  and  the  effect  ofthe  edge  of  the  wing 
from  2).  The  distribution  of  sources  q*  which  we  will  Introduce, 
will  be  necessary  to  create  on  the  lower  side  of  the  wing  a 
distribution  of  pressures  without  accentuated  peaks,  approximately 
llniar,  with  the  exception  ofthe  region  near  the  leading  edge. 


For  this,  for  simplification,  we  will  chose  the  following 
expressions  of  the  contribution  of  sources: 


v'  M 


li 


i x'»»id  i *i oriv'A  i *Hnd  i *.,*))  («><*!<<),  (i3) 

considered  equal  and  of  opposed  sign  for  the  symmetrical 
points  on  axis  0xt  of  the  wing. 


We  obtain.  In  this  manner,  two  "large  wings"  to  form 
a single  one,  having  the  slope  variable  In  such  a way  that 
the  mean  will  be  nought.  For  this  wing,  the  expression  of 
the  axis  of  disturbance  velocity  UlQ  will  be,  slmlltr  with 
(11),  the  following} 


— 


Jo 

cos  h 1 + 


1(1  + >))  |«>«h  1 (1  + By)  ( 1 — Br) 


2B(tj  — x) 


Br>)  (1  + Bs)  | 

B(-t)  x ) 


Jdr, 


2 /•*  ( 

* * ~ \ ’id  + *n»l)[co8h->  1/d  + #»))  (1  — Bj-) 

^ r + 

+ C08h-«I/d  nT-g^Td  . f 

f **,-+*T  jd’,iZ'’ 


(H> 


which  can  b«  put  under  the  form,  after  Integrating, 

;f «•(!•-*•) (cosh' + 

+ 2sin  ' Bl  Bx  I 1 - B*x«  I + *"*«  f B>((*  _ oosU-,  1~  _ 

J 3tt  B»  t B«  _ J.) 


= Qtu  - -*A 
*i  2n  B*  | 


- *•<*•  + x*)  cos  h - Ljbyjf  + 2R*t*  008  h ‘ -i 


.V 


B(t  + x) 


Bx 
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/ 


+ 2(1  - Kl-B*{*)  Bx !l-  B*x*  + 


+ * I At.  [ B'(l'  ~ Tt)  f008  h ■'  + «* h - 1 + 

I 2n  B*  l V.  B(t  _ T)  tty  -f  x)  j 

+ 3J?*<r*coiih-*-i-  + 2(l  - \T-  B*f*)  h - B*x*  j. 

Bx 


+ 


B*(t*  - x*)  cob h_1  ^?  + B»(l»  + **)  oo»h_1 

3nB*[  B(l  — x) 


1 + B*lx  J 
B(t  + i)  *V 


+ (2 »in ' ‘ BIB* »*  + •in'*  Bl-  Bill  - B*/*)Kl  - 


±L. 


Superposltlonlng  the  three  wing  components,  the  resulting 
Imaginary  wing  Is  obtained,  equivalent  from  an  aerodynamical 
point  of  view  with  the  real  wing  for  which  the  axis  of 
disturbance  velocity  has  the  expression 

<ux  = <uu  + <ux,  + nt) 

which  will  be  anti symmetrical  on  the  axis  of  symmetry  Oxj 
and  continue  In  the  origin  0,  We  observe  that  the  velocity 
°11  on  the  higher  side  Is  equal  and  of  opposed  sign  with  that 
on  the  lower  side,  , as  corresponds  to  thin 

wings. 


3.  THE  DETERMINATION  OF  THE  CONSTANTS 

For  the  determination  of  the  constants  q*Q  and  which 
appear  In  expression  (10)  of  Un,  we  will  start  from  the 
following  conditions 
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d» 

da* 

[ dt,  J 

d*  1 

d**1 

l dij  J 

dx  — 0, 


deduced  from  the  theory  of  conioal  motlon(2)f  the  Integration 
being  accomplished  on  a semicircle  d”  of  very  small  radius 
around  a certain  point  y ■ n on  the  wing,  contained  in  the 
interval  (y  » s,  y ■ 1),  (fig.  2). 


We  will  obtain.  In  this  war  the  relations 

, du’ui  i)*  , i\u«,  lt*r*  . 

qi0  ~ a >) ' (7  - h* r”*)3 •* q * a r(  ( i r IF  ^r- ' 

which  stabilize  the  dependence  from  among  the  variation  of 
the  sources  and  the  distribution  of  vertical  velocities. 


Taking  into  consideration  (7),  we  will  start  from  (18) 
and  we  will  put  the  conditions  at  the  limit  in  the  points 
n ■ s and  n ■ i for  the  vertical  velocity  WqjXj* 

wiv-  -taf  {l  - T))n  " HW)-  Kf  - /»’J r*  a,,  Wa) 

l J > IJ* 

- «C  = - ('  ~ TJ><‘  “ #‘V)i  r-  H*  r*  iir„  (Mb) 

HH  ).  T)» 

from  where  we  deduce  the  first  relations  from  among  the 
constants  y*,,  f?  . * 
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—Tgili  — L BI(tia~'  Bl  — Bin-'  Bt)  — 

; — -i-  b* ;* — y y 1 - /?,»,j - 

-ifd- w«-d  -bvH  +ro«h-'-L-co«h-1— U wff 
. 3 l J Bl  Bt\ 


W’.'.V , 
(20a) 


i 

4 2B*l* 


j(l  - B*(*)*/*  4-  [ «*t*(2  - y)  + - 2 yj]  V f - iiSP  + 

+ 3B* /* ( cos  h' 1 — cos  h"'  — ] -f 

l Bl  Bt) 

+ 3«(sin-‘  Bt  - sin*  Bt)  j = »pff  - «*•'. 


(20b) 


On  the  other  part,  the  mean  vertical  velocity  *01x2  or 
the  mean  Incidence  wQ1x2  °f  th«  real  wing,  equal  with  that  of 
the  first  wing  components  as  with  that  of  the  resulting 
Imaginary  wings.  Is  obtained  keeping  In  mind  that  the  two 
large  wings  2)  and  3)  are  compensated  reclprocaly  creating 
a nought  mean  slope;  we  will  be  able  to  write  the  relation 

— ( Wn  t)  <1yi  - tc9Xl, 

( .© 

corresponding  only  to  the  thin  wing  components  1).  Proceeding 
In  this  way  and  taking  Into  consideration  (21)  we  will  have 

y M’ot  l*  = ~ •C  *’  + j m’o.  vH  (22) 

which,  after  accomplishing  the  calculations,  becomes 
{»  [**•»-•  Bl  - star*  Bt  + Bl[l  - B*Py*  - B#(i  _ 4. 

J + 2 [ Bl(  1 - B* l*)**  - B,(i  - /?«*i)»'*|  4 

M 

f ‘ + - <!  - «*••)•*)}  - wiV-  «r„.  (23) 

t 


15. 


Ve  will  determine  the  constant  A21,  which  appears  In  the 
expression  (9)  of  the  axis  of  disturbance  velocity,  considering 
the  variation  of  vertical  velocities  w^j  from  a point  on  the 
wing  to  the  one  of  the  nought  vertical  velocity  (for  example  on 
the  Mach  cone).  In  order  to  avoid  the  difficult  calculations 
which  appear,  we  will  consider  that  the  sources  distributed 
on  the  line  In  the  Interval  (s,l)  are  concentrated  In  y » n' 
of  Intensities  Q20  and  Q2I  in  such  a way  that  we  will  have 

= ~ ( 1 ~ y ) ’ “ 4"  9*‘ 1 ( 1 “ 7)  »•) 

Prodeedlng  In  this  way,  we  will  write  the  relation 

rul  Mvr=rHrzk*<it'u  , 

* J TT  dx  = »#> » IIS' 

•VrlpA  JL 

where  Is  the  axis  of  disturbance  velocity  for  the  simplified 
dase  of  the  source  concentrated  In  y * s*,  given  by  the 
expression  

- -r  *;  - + t « 0. + v«  ■ ')  h ■ 1 - 


2 . ..l/ir+V) 

(<^M  - V11  •*)  <*>8  h —£[77 


V)(t  + J) 
(«'  + »•) 


Accomplishing  the  Integral  (25)  on  the  axis  of  the  ordinates 
m 0,  x ■ 1^)  between  the  limits  0 andoc,  results 


(2  — B*l*)E(k)  - H*l*K[k) 
1 i»(i  — B*l*) 


+ JLq:J  1 _ 1 (1  - B'»')E(k) 

4«  l i J«'l L 
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- «»<■(£(*)  _ ll(p,jt»j+  i 

x _H_  / _t  /.  »'  vj  b*  Yi» i'i 

**,ll  ~T)  T~^l's ■-■'I 

in  which  K(  k) , E(k),  II(p,k)  represent  the  complete  ellptlcal 
integrals  of  the  first,  second  and  third  Instances  respectively, 
having  the  module  k and  the  parameter  p given  by  the  relations 

..  .....  . K1  I n - K(k)K(<v0,n  + 

1Up.*>  MM  i J^yjrrv*  [ -2 

1 i:Vl  ' 

(K  (A)  - KlA»m..*  > • 


I | I - JPP,  p />‘J«'a  I - NHI  I -*)*  - ». 
A’  /»/,  9o  »»»  '—• 


l -'«»>) 


Furthermore,  due  to  the  separation  of  flow  at  the  edges, 
we  will  Impose  the  condition  that  the  velocity  be  finite  at 


the  subsonic  edges,  canceling  the  constant  A21* 


A21  - °* 


(29) 


Eliminating  *^f)  and  w*J]  between  the  equations  (20a), 

* . * 

(20b),  (23)  and  (27),  we  will  obtain  the  constants  q20  d2l* 


V;'’  / i i i **’•»»  v?i  *— — ((■„„  <ao) 

31'  30  ^ao'^HI  / i/ — /w>/|, 

In  whloh  we  made  the  notations 

1 + i-jyi  _ _ 2h  - B*t*  - A JU(al in  • IN  - 

sin-1  H»)  - -1-  P (K  1 - R*P  - — V 1 - #»»«»)  * ( 

- — [ (1  - fPP)w  - (1  - «••*)"*]  + cosh  • — -M«h  » — J, 
3 Bl  Ht  | 

17* 


U ' (1 " K'P)‘ " + - f ) + £(‘ -af)]inrw?  + 

+ 3BlJ»in-'  Bl  - ,i„-i  B.  + Bl^foah  ' J-  - conh-'i-J], 


(31b) 


g«r(3  [8in~,W  - Bin'1  B*  + B/( 1 - Btyy*- 
+ 2 [Bl(l  - fit  {.)./,  _ B#(1  _ B,,l)>/t]  + 


B»(I  - /?*#«)'  *]  + 


+ 6B/[<1  ~ ~ U - B* 


**•'{  l ) I - fcv* 


»*)•'*]  J 


(I  - b*#'*) B(fc)  -Bv*<ir<*,  _ n<P,  *» 

(31  <*) 


- - i h f 11  - *“■>* ! + [ "■  *■  (*  - j)  ■ + > (*  - » t)]^-  . 

+ 3 B*  l*  ( cos  h ~ 1 — — cos  li _ * — ) 4. 
l Bl  Hh) 

1 Bl  — sin-1  Bs)  | 4- 


+ 3 Bl  (sin" 


B»(*^l  - ' ' 

^(i  - 7)'— l<a  - nic u - K(H).  mil 


Th«  constant*  kjQ#  k20»  *n  an<*  ^21*  whlch  appear  In  the 
expression  of  Ulc  given  by  (14),  Is  determined  taking  Into 
consideration  the  role  of  the  third  wing  component  which  will 
have  mean  slope  -w0lx2|  similar  with  (21)  we  will  write 
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J Jo 


from  where  we  deduce  the  relations 

• I 

(«.,  - ri -•  i„l  n < i * *l0n>(i  - n-)** **n» 

Jo 

(«•„  - - *»  ('  <M  *„r()(  1 - «•>>»)“  *‘n. 

#**  J * 


l -l.'*) 

(WO 


which  ends  the  form 


M'iV  — "’tfl 


-***-  Jl((i  _ /;»/*)“  _ i)  (i  - /W- 

- ;{ 


•4V  - •<', 


1 i I,  — , „ . sin  1 HI 

--(I  — H*l*)**  - l I HU*  f 

1J  H \ HI 


k"  U »/  •_•<»  ft*!*)**  ) all  r-  /<>/«  + — 

1W  \ h l W 


— ^ |(t  - HU*)---  ijl- 

Considering  for  kt0  and  kn  suitable  ralues 


l|tl  ~ 1|1 


1 


i;u»> 


we  will  obtain  from  equations  (23),  (34a)  and  (34b)  the 
constants  k2o  *nd  k21» 

«*  j f 

«fwj  3 Jfr(i  _**•).,,  + 

V; 

4-a  [Si(i  - - />« (i  _ »•*•)*»]  ^ 

*•*)••  ]Jx 
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S’ 


*■  L 10  12 


B»t*  + 


•in  ~tBl\  4 1-* 


HI 


(36a)  > 


^ ~ 9"  | 3|8in-*Bf  - + Bf(l  _ BM*)1'1  _ 

, + 2|b/(1  — B*/*)w  - B#(l  _ B* «*)*'*j  + 

P 

+ ~un  [ (1  “ “ (1  ~ ^’**)5/1]J  x 

{t  B#p[2(1  -«*‘,)w+3|Kr=~Bni  + -ip^  BiJJ  + 


(36  b)1 


4.  THE  DISTRIBUTION  OF  PRESSURE  AND  AERODYNAMIC 
CHARACTERISTICS 


We  have  shown  that  the  axis  of  disturbance  velocity  on 
the  real  wing  results  through  the  superpositioning  of  the 
three  imaginary  wing  components,  obtaining  formula  (16).  In 
this  way,  the  coefficient  of  pressure  will  be  obtained 
considering  the  axis  velocity  given  by  (16),  (fig.3): 


c = _ 2 = 

vm 


2Re 

Um 


(37). 


For  the  calculation  of  the  coefficient  of  lift  of  the 
straight  line  of  the  wing  we  will  make  the  observation  that 
the  wings  of  symmetrical  thickness  do  not  give  lift,  so  we  will 
only  consider  that  given  by  "the  lift  wing": 
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1 


which,  after  accomplishing  the  calculations,  terminates  the  form 


Fl«  * 


tvr 


It 

I 


For  the  definition  of  the  parameter  j»,  which  enters  the 
expressions  above  and  which  determines  the  limits  of  the 
distribution  of  sources,  we  will  remark  first  that  the 
position  of  maximum  distribution  of  pressures  colnolde  with 
the  abscissa  y ■ o of  the  oenter  of  the  vortex  nucleus,  as  Is 

22, 


Lk 


4 


ascertained  from  experience.  However,  making  the  calculation 
on  the  base  of  distribution  of  sources  (6),  It  Is  ascertained 
that  the  peak  of  the  depressions  on  the  higher  side  of  the 
wing  falls  approximately  at  half  the  distance  between  the 
center  of  gravity  of  the  Intensity  of  the  sources  and  the  point 


of  abscissa  s,  which  we 


follow  to  determine  It.  In  this 


way  we  can  deduce  the  relation  between  s and  c 


«•  K f (I  — S) 

ti 


1 0.-‘ 

I 


For  the  definition  of  the  positions  of  the  vortex  nuclei, 
we  will  observe  that  at  small  torsion  ofl  the  wing,  this  can  be 
considered  approximately  plane  and  parallel  with  the  direction 
of  the  undisturbed  flow  U*o  , so  that  the  vortexes  falling  off 
at  the  leading  edges.  In  the  form  of  a horn  would  be  considered 


as  attached  to  thtfvsurface  of  the  wing*  This  work  would  lead, 
however  small  the  Intensity  of  the  vortexes  would  be,  to  very 
great  local  velocities.  Incompatible  with  the  real  effects  of 
the  separation  of  flow  at  the  edges.  In  order  to  avoid  this 
work  we  must  allow  that  at  very  small  torsion  of  the  wing, 
therefore,  for  very  small  values  of  the  parameter  <X01  ( 2 ) or, 
more  exactly,  when  0,  the  position  of  the  vortex  nucleus 

will  be  c ■ 1. 


For  greater  torsion  of  the  wing  we  will  allow  that  the 
nature  of  the  vortexes  which  start  on  the  leading  edges,  evolve 


proportionately  with  the  Incidence  ol  ■O^01x2  as  BS  with 

the  span  of  the  wing  (as  with  the  #tng  with  constant  Incidence 
(1)),  therefore  with  the  square  of  the  span  of  the  wing,  and 
winds  Itself  on  an  axis  representing  the  line  of  the  cente^ of 
gravity,  corresponding  to  the  abscissa 


c 

1 


We  will  remark  however  that  the  Intensity  of  the  vortex 
nucleus  Is  something  smaller  than  the  total  vortex  generation 
Intensity,  the  rest  being  In  the  layer  which  forms  the  surface 
of  the  horn  along  the  leading  edge. 


From  that,  orienting  ourselves  according  to  the  experiences 
accomplished  on  the  delta  wing  with  constant  Incidence,  the 
minimum  position  of  the  vortex  center,  corresponding  to  torsion 
or  greater  parameter  w0j  ■ -(Xoiu*>*  can  b®  considered  to  contain 
between 

<•  as  o,ti  / -m>,7  /,  (a) 

that  which  would  correspond  with  the  center  of  gravity  of  only 
a part  of  the  total  Intensity  of  the  vortexes,  representing 
Just  the  real  Intensity  of  the  remaining  nucleus  (about  80<) . 

What  makes  this  minimum  position  to  be  the  weight  of 

Unter_laEJ 

stability  Is  the  fact  that  at  gre^t  incldenceSvvortex  nuclei 
appear,  of  opposed  sign  with  the  principal,  situated  between 


this  and  the  leading  edge 


Between  these  limits,  taking  into  consideration  the 
experimental  results  obtained  by  various  authors  on  tbe  plane 
delta  wing  with  constant  incidence  we  will  allow  the  following 
approximate  formula  of  variation  with  the  incidence  for  the 


position  of  the  vortex  nucleus: 


/ 1 + 1,1  (a)1'- 

where  oc  represents  the  incidence  of  the  wing  between  a suitable 
point,  representing  the  center  of  gravity  of  the  aerodynamic 
effects,  namely 


;i  _ iL  j /=  — / 


(Hi) 


We  will  remark,  that  the  whole  reasoning  which  has  led 
to  the  stability  formula  (45)  can  apply  for  each  section  of 
the  wing,  obtaining  in  this  way  a curved  line  for  the  position 
of  the  vortexes  nuclei.  However  the  deviations  from  a straight 
line  is  ascertained  toward  the  peak  of  the  wing,  where  the 
aerodynamic  contributions  are  very  small,  and  behind  the  wing, 
where  the  incidences  grow,  the  vortexes  sit  rectlllnlar.  From 
this,  we  can  allow  a straight  line  for  the  axis  on  which  the 
vortex  sheet  winds  in  the  form  of  a horn,  corresponding  to 
the  section  given  by  (46). 


This  hypotesls  permits  us  to  apply  the  methods  of  conical 
flow  of  the  higher  order. 


25. 


5.  THK  CASE  OF  CONCENTRATED  SOURCES 


Considering,  for  simplification,  a sudden  variation  of 
vertical  velocity  in  the  point  a/ • o,  corresponding  to  the 
abscissa  of  the  vortex  nucleus,  we  will  obtain  the  case  of 
an  "edge"  of  separation  of  vertical  velocities  w^cxj  and 

cxj  (2). 


The  sudden  drop  of  vertical  velocity  or  of  incidence  is 

realized  analytically  with  the  help  of  logarithmic  singularities 

concentrated  in  the  points  s'-  a c on  the  track  of  the  wing 

in  the  auxiliary  plane  x (4),  having  the  intensities  Q „ and  Q,.  . 

20  c 1 


In  this  case,  the  axises  of  disturbance  velocities  of  the 
three  wing  components  have  the  following  expressions 


•i*. 


IV-.1I  V,|  .I1)  VON  ll 


J'  ' .+  “ IV„,  i •]/»'*  <>«-•*) 

i*  * r ^ j) 

f'dv 

t ; ( i v»  h ■ ; i v..  ^ u - 

for  the  "thin  lift  wing"; 

~ (V-,,  I V»|J)  «•»•*  li  'I  " 1 ' 1 ,h) 

•'«  71  l •-'/{(«•  .,) 

(V»,,  V,l  •*')  tl'N  ll  'I  *'  t /•*<•)  (|  I IU)  ~ 

* i - /•(<•  ( .r) 

“ ( I V*>  *'»•*  l»  ' If  ' li>rt  , v„.<  .„sl,  'I  ' . 

” v r /<*(«•»  .r»)  | ,4,1  * i 
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lit) 


• i 


for  the  "large  wing" , and 

•Ml  ^0 

for  the  "compensating  wing"  for  slope,  where  (Uu  has  the 
expression  given  by  (15).  Given  that  we  will  remark  that  the 
sign  of  the  axis  velocity  Is  different  on  the  two  sides 
of  the  wing.  The  expression  of  the  total  axis  velocity  will 
be  given  by  a formula  similar  with  (16),  which  will  permit  us 
to  determine  the  distribution  of  pressures. 


Starting  from  (17*)  end  (l?b),  we  will  find  In  the  case 
of  concentrated  souroes  for  Q2o  ®nd  ^21*  the  folloing  expressions! 

w“ - **»•*»  ^ 

and  from  (25),  In  which  we  Introduced  <M\,  given  by  (4?),  we 

will  obtain 

. (i!  — H,l*)K(k)  — 

**t\  .....  

i*  (l  — 11*1*) 


ft  W« 


y?  c* 

( I - It*  <•*) 


[I  H*c*)  L (k)  - Htl*{K(k)  - II  (j>,  till 


•'  It*]/ 1*  -<•*  \ 

+ - j-it*,»  l(‘  ~ ll‘ ,,,) 11  <?• k)  - K (*>l  <'• 


In  which  the  complete  ellptlcal  Integrals,  In  this  formulm, 
have  the  module  and  the  parameter  P given  by  the  relations 
(28a)  and  (28b),  In  which  will  be  put  % - o. 


Using  the  relation  (21),  we  will  obtain 

(«c  - - (»*•„,  - 
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„ V 


which,  together  with  the  condition  that  the  velocity  be 
finite  at  the  edges 


*2i  ■ 0» 


(53) 


and  with  the  equations  (50),  and  (51)  fora  the  system  of 
equations  from  which  we  deduce  the  constants. 


.• 


«V> 


In  this  manner,  we  will  obtain  respectively 

£(1)  - B*l-  [K  (k)  - (l  - H* <•*)  11 1 p.  1)1 


A’<*)  - B'  P [K  (k)-(l- P*e»)  n(p,l)  1-^(1  ('*-  ll(i-u*f*) 


E(k)  -£V*)Il(p, *)]- 


it(l  x 

= (..  £)l!  'I" 


£ (1) - B‘  i*[  £ (1)  -( 1 - Bt <•*)!!( p,l) j - * ( 1 - £'0)|  (--!)(!-  **,-*) 


<,>»  = ~ 


(.VI  a) 


(54  l>) 


- ~)l  1 -£*<•* {£(!)- (JCl*)  - (1  - B*c*)  II  (P,l)  | - i d - B*c*)  [/  _j)  (1  _ j*.^ 


(54  c) 


£2  r3 


2 |i-^jVi-n««»j£(*)-£tt*[i:(t)-d-B,c»)n(p,fc)]-^-d -»•«•)  / (£_i)d  - i»*p*)|(»4d) 


The  constants  k2o»  k deduced  In  (34a),  (34b),  in  which 
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w^  Is  substituted  with  the  expression  given  by  (54b). 


For  the  determination  of  the  coefficient  of  lift  of  half 
of  the  wing,  we  will  start  from  (38)  In  whloh  the  expression 
ufji  ■ Re^*'  will  be  Introduced  from  (47)  and  we  obtain: 


28. 


r 


‘ I + 7*0<"b  ‘^)) “ 

illr.j  - H*  .-*)«»  It  * - - i*W*  1/  l - 
~ a*  ('*  / I r r 


c* 


(W 


The  coefficient  of  the  moment  of  roll  le  calculated 
starting  from  the  relation  (40) 


lie. 


’,c-b-+ 


(<V’  j x £1  h _ /;*  (i*  + tv*)  j . 
,'i  I * l 


<Mi> 


In  which  the  constant  Q2q  has  the  expression  (54c),  given  above. 


6.  FINAL  REMARKS 


The  results  obtained  above.  In  the  case  of  toreloned 
wings,  can  be  used  with  the  plane  delta  wing  In  rotating 
motion  around  the  axis  of  symmetry  Qxj. 

In  this  oase  however,  due  to  permanent  rotation  of  roll. 
It  will  be  necessary  therefore  to  make  other  considerations. 

On  the  other  part,  we  will  remark  that  the  position  of 
vortex  nuclei  will  be  affected  by  a centrifugal  effect  , which 
pushes  them  towards  the  leading  edges  removing  them,  in  the 
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case  of  great  velocities  of  rotation,  to  greater  distances, 
even  over  O.75  1, 
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